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Abstract 

Experiments designed to measure neutrino oscillations also provide major opportunities for dis- 
covering very weakly coupled states. In order to produce neutrinos, experiments such as LSND 
collide thousands of Coulombs of protons into fixed targets, while MINOS and MiniBooNE also 
focus and then dump beams of muons. The neutrino detectors beyond these beam dumps are 
therefore an excellent arena in which to look for long-lived pseudoscalars or for vector bosons that 
kinetically mix with the photon. We show that these experiments have significant sensitivity be- 
yond previous beam dumps, and are able to partially close the gap between laboratory experiments 
and supernovae constraints on pseudoscalars. Future upgrades to the NuMI beamline and Project 
X will lead to even greater opportunities for discovery. We also discuss thin target experiments 
with muon beams, such as those available in COMPASS, and show that they constitute a powerful 
probe for leptophilic PNGBs. 
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I. INTRODUCTION 

In the last decade, several experiments have explored neutrino masses and mixings, but 
these high-luminosity laboratories are also sensitive to rare production of new metastable 
particles. Neutrino beams such as the LAMPF Neutrino Source at Los Alamos and NuMI 
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and BooNE at Fermilab are produced through two basic stages. First, a high-intensity proton 
beam impinges on a target and produces a large number of pions (and other hadrons), which 
decay to muons and neutrinos. The muons are stopped in a thick layer of rock, while the 
neutrinos travel unimpeded through the rock to the detector. In fact, short-baseline neutrino 
detectors are situated behind the most intense proton and muon beam-dumps to date. Thus 
they are ideally configured to search for long-lived particles produced by rare proton-nucleus 
or muon-nucleus interactions. 

Two classes of new physics scenarios naturally give rise to light, feebly coupled particles 
of this type. An approximate symmetry broken at a high mass scale F naturally gives rise 
to light pseudoscalars — pseudo-Nambu-Goldstone bosons (PNGBs, or "generic axions") 
- with couplings of order mx/F to Standard Model matter X. Alternately, a new "dark" 
U(l) gauge boson can naturally have small kinetic mixing e with the photon, giving rise 
to suppressed interactions with all electrically charged matter pQ. Either spin-1 or spin-0 
bosons can be radiated in energetic-particle interactions with matter, with a very small rate 
proportional to the square of their weak coupling. The luminosities achieved in fixed-target 
experiments are such that thousands of these particles could be produced, and they offer a 
new window into weakly coupled sectors. Once produced, the lightest particles in a hidden 
sector could decay only through their weak couplings to Standard Model particles, and so 
would be quite long-lived and weakly interacting. While ordinary products of the collision 
are stopped in the shielding upstream of a neutrino detector, exotics could penetrate the 
shielding and decay within the detector, yielding a distinctive signal. 

This note summarizes the several classes of exotic particles that naturally give rise to 
observable late-decay signals, their experimental signatures and typical kinematics. In any 
given model, the production cross section and lifetimes of these exotica are both determined 
by a single small coupling parameter and by the mass of the produced particle, so in par- 
ticular, we present specific estimates for the sensitivity achievable with late-decay searches 
in MINOS /MINERvA, MiniBooNE, and LSND. 

PNGB's coupled to hadrons were searched for extensively in both proton and electron 
beam-dump experiments in the 1980s, most notably in CHARM [2J at CERN, E774 [3] at 
Fermilab, and the SLAC experiments E137 [1] and E141 [5]. Many of these limits have 
recently been re-interpreted ([SHE]) in the context of kinetically mixed gauge bosons and 
the associated scalar bosons that give them mass through the Higgs mechanism. Kinetically 
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mixed gauge bosons have been a subject of considerable recent interest [9H2T] and discussions 
of other collider, accelerator, and direct and indirect astrophysical probes for them can be 
found in e.g. In particular, [7J discussed the sensitivity of neutrino experiments to 

hadronic production in the case of kinetic mixing and the potential importance of existing 
LSND data as a constraint on these models. Additional constraints on both classes of models 
from supernovaes, rare decays, and radiative corrections have also been extensively discussed 

PUS]. 

Our aim in this note is to present a more complete analysis of the sensitivity of past, 
present, and future neutrino experiments to new weakly-coupled physics, with a particlar 
focus on PNGB models. We hope such a unified summary will facilitate new analyses of 
neutrino-detector data to discover or constrain new weakly-coupled particles. We also dis- 
cuss the potential reach for experiments with muons beams that strike a fixed thin target. 
In ^TTj we review constraints on generic pseudoscalars from other arenas, in particular su- 
pernova data, rare meson decays, and the anomalous muon magnetic moment. The range 



explored by neutrino experiments is complementary to all of these. In { III we consider the 
implications of existing LSND analyses for both PNGB's and dark gauge bosons. Due to 
the large number of protons dumped in LSND, we find that these analyses provide stronger 
constraints on PNGB's than the CHARM experiment. Our results for dark gauge bosons 
are consistent with [7], but we have tried to clarify the experimental sensitivity. In £IV we 



consider the sensitivity to PNGBs that could be achieved by analyses using modern neu- 
trino beamlines, such as BooNE and NuMI, and their near detectors (the complementary 
analysis for kinetically mixed gauge bosons was presented in [7j). We consider both the 
standard production mode in proton-nucleus collisions and the production of very forward 
PNGBs off the stopping muons. The second mode is enhanced by the magnetic focus- 
ing of pions, so experiments using focused neutrino beams are uniquely sensitive to purely 
leptophilic PNGBs that do not couple to hadrons. For ordinary PNGBs coupled to both 
quarks and leptons, searches in MINOS /MINERvA and MiniBooNE would have sensitiv- 
ity comparable to, or perhaps slightly better than, the CHARM beam-dump limit [2\. A 
future "Project X" could significantly extend this reach into new territory. For leptophilic 
PNGBs, MINOS/MINERvA has slightly better sensitivity than the constraint from E137. 
In §[V] we discuss muon fixed-target experiments using thin targets. Such an experiment 
could be possible at the COMPASS experiment at CERN. We find that a COMPASS-like 
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setup can probe new territory in the parameter space of leptophilic PNGBs, closing a gap 
between muon g — 2 limits and those from neutrino experiments. We conclude in SVI, and 
an appendix discusses the details of PNGB production off muon beams. 



A. Models of Weakly Coupled Light Exotics 

1. Pseudo-Goldstone Bosons 

Light pseudoscalars can arise as pseudo-Goldstone bosons in a large variety of well- 
motivated theories, such as multiple higgs doublet models, theories with an R-axion [45] 
(from spontaneous and explicit R-symmetry breaking in a supersymmetric theory), axion 
models [46J, the Next-to-Minimal Supersymmetric Standard model [47J (NMSSM), and re- 
cent dark matter models with a scalar portal to the dark sector [IS]. The most important 
point is that these particles are naturally light if there is an approximate shift symmetry, 
and that their interactions are proportional to the inverse of some symmetry breaking scale 
F. Using fermion equations of motion, the derivative coupling of a PNGB a to a fermion 
bilinear turns into the operator (in Weyl notation) 



which is the coupling we assume for the leptons and/or quarks, Xi °f t ne Standard Model. 
A phenomenologically interesting sub-class of PNGB models are those where the PNGB is 
leptophilic, i.e. it couples preferentially (or only) to leptons; this scenario could arise if the 
lepton sector has its own higgs mechanism separate from that of the quark sector. Leptophilic 
PNGB are essentially unconstrained by searches for rare meson decays and proton fixed 
target experiments, so experiments that have muon beams, such as MINOS/MINERvA, 
MiniBooNE, and COMPASS, can easily be the most sensitive probes of these particles. 

The coupling of a PNGB with mass m a to a fermion with mass rri£ induces a PGNB 
partial width 



and the total width is well approximated by T e + V ^ for m a < 400 MeV (for larger masses, 
hadronic decays can also become important but we use the leptonic widths for masses up 




(1) 




(2) 
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FIG. 1: Left: The rest-frame lifetime of a pseudo-Nambu-Goldstone boson (PNGB) as a function 
of its mass m a and decay constant F. Right: The lifetime of a dark photon A' as a function 
of its mass m^' and e, the strength of its mixing with the Standard Model hypercharge gauge 
boson. In both plots, the black lines correspond to different decay lengths (cr): 10 /im (solid), 1 
cm (dot-dashed), 1 m (dashed), and 100 m (dotted). In the blue, purple, red, green, and white 
shaded regions the decays are prompt (< 10 /mi), displaced with < 1 cm, displaced with > 1 
cm, "invisible" with > 100 cm, or "invisible" with > 100m, respectively. In fixed-target or beam 
dump experiments the particles typically get a large boost that increases their decay length by 
-^beam/ m ass. The feature in the left plot at 2m M occurs since PNGB's coupling to a Standard 
Model particle is proportional to that particle's mass, and at this point decays to two muons are 
allowed. The dip in the right plot near 0.7 GeV is due to the p-resonance. The lifetime for both 
the PNGB and the A' is calculated assuming decays directly into Standard Model particles. 

to 1 GeV). Thus, for example, proper lifetimes of 1 mm are obtained with F « 70 GeV 
(m a = 100 MeV) or F w 20 TeV (m a = 300 MeV). Fig. [j] (left plot) shows the decay length 
(cr) of a PNGB as a function of its mass m a and decay constant F. Note that the decay 
length is very different above and below the muon threshold, due to the much stronger 
coupling to muons as compared to electrons. We see that for F < 10 2 GeV, they decay 
promptly and colliders should be able to set the best constraints. For larger F, collider 
searches that look for displaced vertices or missing energy can still set limits, but searches 
in beam dump experiments (with a large shield) become relevant. 

We ignore decays of PNGBs to two photons, since this is always subdominant in the mass 
range we consider in this paper. 
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2. Kinetically Mixed Gauge Bosons 



Another class of light particles that has received significant interest in the last few years 
is a light, weakly coupled "sequestered sector" lying alongside the Standard Model. New 
states at the MeV-GeV scale are not in conflict with data, because the gauge and global 
symmetries of the Standard Model greatly restrict the couplings of ordinary matter to these 
states. It is, however, quite natural for the few allowed interactions to be suppressed by loop 
factors. 

In particular, we consider a new MeV-GeV scale Abelian gauge boson A' coupled to 
electrically charged Standard Model particles ip 

5C = eeA'^ij). (3) 

Such a coupling can generically originate from the kinetic mixing between the field strengths 
of the Standard Model hypercharge and a hidden sector U(l) PQ, 

SC = |W (4) 

where F' = d^A' u — d u A' is the field strength of the A' gauge boson, and similarly Fy V is 
the hypercharge field strength. This mixing, assuming the hidden U(l) is broken so that 
the A' is massive, is equivalent in low-energy interactions to assigning a hidden charge eeqi 
to Standard Model particles of electromagnetic charge where e = ey/ (cos 9\y) and 9w is 
the Weinberg mixing ang le (Eq. @ if q f = 1) . Such a mixing can be generated in many 
ways, e.g. through loops of new heavy particles that couple to both the A 1 and Standard 
Model hypercharge. We refer the reader to e.g. [27] for more detailed discussions. 

In Fig. [I] (right), we show the decay length of a vector boson A' as a function of its mass 
rriA' and the parameter e which sets the strength of its kinetic mixing with hypercharge. 
Assuming the A' decays into Standard Model particles rather than exotics, its lifetime is 

3 _ 0.8cm / E Y 10 ~ 4 ^Y 100MeV 

where we have neglected phase-space corrections, and N e g counts the number of available 
decay products. If the A' mixes kinetically with the photon, then A^ fj = 1 for tua' < 2m M 
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TABLE I: Summary of Signal Properties for Light Exotics considered in this paper. The PNGB 
decay to two photons is never important for the mass range we consider, and can be ignored. More 
complicated signals are possible if the PNGB or A' can decay to other hidden sector particles before 
decaying to Standard Model particles, but we will not consider this possibility. 



when only A' — > e + e decays are possible, and 2 + R(rriA') for tua' > 2m M , where R = 

o-(e+e--> liadronS;g=TO A /) rjqi 
o-(e+e--^+/i-; E=m A ,) 1=^1 " 

We summarize the possible production and decay mechanisms for both PNGBs and A n s 
in Table |lj In this paper, we only consider the case for which the PNGB or A' decays directly 
to Standard Model matter. More complicated signals are possible if they can decay to other 
hidden sector particles before decaying to Standard Model particles. 

II. CONSTRAINTS ON PNGBS FROM RARE MESON DECAYS AND SUPER- 
NOVAS 

A. Limits from Kaons and B-Meson Decays 

In this section and in Table [TT| we briefly summarize constraints from meson decays 
on PNGBs with sub-GeV masses. These constraints are controlled by two factors: the 
partial width for the rare meson decays into a PNGB, and the fraction of PNGB's that 
decay promptly enough to be included in the data samples (or, in the case of invisible-decay 
searches, the fraction that decay outside the detector). The combination of these searches is 
sensitive to PNGBs with F < 10— 100 TeV over a wide range of masses. We call these limits 
"exclusions" below in the interest of brevity, but it should be emphasized that several of 



S 



Decay mode 


B.R. limit 


Mass range 


Decay region 


R+ — i. K+ -1- •inn) r^Til 
D — 7 1\ ~r LllV. lOUl 

K + -)■ 7T+ + m«. [5J] 

-> 7T+ + X X -> e+e" [52! 
5+ -> if+e+e" [53] 
dT/dm 2 (B + -»■ K+£+£-) [M] 


0.5 - 10 x 10~ 10 
1.5 - 5 x 10~ 6 
6.7 x 10~ 7 
1.2 x 10~ 7 


U v. //t a <^ ^liLfM 

< m a < 2my 
10 MeV < m a < 120 MeV 
30 MeV <m a < 2mJ i * ) 
140 < m a < 1440 MeV 


^> l.UOItl 

> 1.3m 
< 5cm^ 
< lOO/W^ 
< .5cm 



TABLE II: Summary of constraints on PNGBs from various experiments (B.R. = Branching 
Ratio). For decay of meson X, a characteristic transverse energy mj/2 is assumed in computing 
the fraction of decays either in or outside a cylinder whose radius is given in the "decay region" 
column. 

M; Searches extend beyond 2m M but only imply relevant limits on PNGB models below 2m „. 
ft); Length scales are guesses only, but overall exclusion is insensitive to cutoff because it overlaps 
with reach of invisible searches. 



the results have been re-interpreted by non-experts (the authors) in a context very different 
from the original experimental design, from which significant inaccuracies in the boundaries 
could have resulted. We focus here on the most constraining searches; a more exhaustive 
set of limits are considered in [43j, but our treatment of the decay regions and the PNGB 
lifetime differ from theirs. 

Following [32]) we consider two-higgs-doublet models as a generic framework for a new 
PNGB a that couples to the Standard Model (whereas [33] assumes NMSSM couplings). 
In this framework, rare meson decays to a lighter meson and a are mediated primarily by 
a top-quark loop, which receives contributions at all scales up to the highest energy at 
which the PNGB couples to the top quark (typically either the F scale or the electroweak 
symmetry-breaking scale). Therefore, the meson decay rates depend on the detailed UV 
dynamics through which the a particle couples to Standard Model fields, as discussed in 
j. We adopt the approximate formula 



T{B + -> K + a) 



G F \V tb V ts \ _ m 4 m | (^^j (kinematic)[f (m a ) 2 } 



V22 



(6) 



and likewise for K + — > 7i + a, but with the product of CKM matrices V t * b V ts instead given 
by V t * s V t d- Here Gp is the Fermi coupling constant, v = 174 GeV is the electroweak Higgs 
vacuum expectation value, m t is the top mass, and m# the 5-meson mass. This formula 
is obtained from (A.l) of [32] by assuming (3 = 45° in the two-higgs-doublet model, replac- 
ing sin^ — > v/(2F), and setting the model-dependent combination \X\ + X 2 \ = 1. This 
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combination depends on the physical Higgs mass but is typically in the range 1-10 so that 
our choice is conservative (although note that there are special choices of the Higgs mass 
where the sum of X's crosses zero, in which case limits from meson decay are weaker). The 
kinematic factor is typically in the range 0.5-1 and is defined in (A. 3) of [12], while the 
form factor /q w 0.33 [55]. The formula above with our parameter choices yields branching 
fractions about a factor of 10 smaller than those in [15] . 

The second factor that comes into play is the decay length of the PNGB. For visible 
meson decay searches, the PNGB must decay within a small distance (typically of order 
mm-cm) from the meson decay from which it originates, so that all of the tracks produced 
by the meson reconstruct a single vertex. Likewise, invisible-decay searches require that the 
PNGB decay outside the detector volume, with a typical size of about a meter. We have 
attempted to take these effects into account more precisely than [13], which uses a uniform 
length scale of 10 m for all limits. 

The PNGB partial width to £ + £~ is given in Eq. ^ and the lifetime is shown in Fig. [Tj 
Naively, one would expect that PNGB's with characteristic lifetimes between 1 mm and 1 
m might be poorly constrained by both visible-decay and invisible-decay searches. However, 
the fraction of events in which the PNGB decays within an atypically short distance from 
the meson vertex can still set significant constraints; moreover, the limits from experiments 
with different energy scales (mjc and m#) overlap to fill this intermediate-lifetime region. 

a. Invisible decay searches for B + — > K + vv or K + — > n + X, where X decays invisibly, 
peak in a specified mass range. Below the muon threshold, these are sensitive to the highest 
F's probed in accelerator experiments. 

• B + — > K + + inv.: A search at Belle [50J set an upper limit of 1.4 x 10~ 5 on the 
rate of B + — > K + vv. This search imposes an upper limit on the energy deposited in 
the electromagnetic calorimeter (ECL), which extends to 1.65 m from the beamline. 
Therefore, for this limit we consider only a decays outside this radius, assuming a typ- 
ical transverse momentum of mB/m a . A tighter limit could likely be set by searching 
for a narrow invisible resonance in the B + decays. 

• K + — 7- 7r + + inv.: The search for K + — > 7t + vv at Brookhaven E787 [51] also set an 
explicit limit on the decay K + — > n + X where X is invisible. The branching fraction 
for this mode must be below about 5 x 10~ n for X masses below about 100 MeV, and 
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about 10~ 9 for X masses between 150 and 250 MeV. No limit is set between 100 to 
150 MeV, where several K + — > n + uu candidates were seen. The transverse size of the 
detector is roughly 1.3 m, and the kaons are produced at rest so that m K /2 is a good 
approximation to the energy of the outgoing PNGB. 

Of these two, the K + decay limit extends to higher F but the B + decay search is able to 
probe lower F because the PNGBs produced in B decays are more boosted and therefore 
longer-lived. There are additional constraints from CLEO and BaBar searches for T(ls) 
and T(3s) decays to 7a [56H58] . but this region is largely contained in the two identified 
above. 

b. Visible decays Again, we focus here only on the most powerful visible-decay searches: 

• K + —7- n + X, X — » e + e~: This resonant decay mode was searched for in the 
experiment at KEK [52], which excluded kaon branching fractions to below about 
1.5 x 10~ 6 for PNGB masses between 10 and 80 MeV, and below about 5 x 10~ 6 
for PNGB masses from 80 to 120 MeV. The experiment also set a limit above 140 
MeV, but it is not as constraining as the -B + -decay limits discussed below, so we do 
not include the constraint in this region for this mode. used stopped kaons, so 
we assume an initial energy of m^/2. No vertex requirement is explicitly mentioned 
in the analysis, but common vertex requirements of order a few cm are frequently 
imposed in spectrometer analyses typical in these experiments; in any case, at much 
larger distances, the mass resolution would likely be degraded. We have used an ad 
hoc but conservative estimate of 1 cm to produce the limit in Fig. [2] but even if the 
vertex requirement is much tighter (as tight as 0.5 mm) the excluded region would 
overlap that of the B + — > K + + inv. search. 

• B + K + £ + £~: BaBar [53] and Belle [5U [59] have both measured the rate of the rare 
decay B + — > K + £ + £~ for £ = e, /i. The observed branching fractions, (3-6) xl0~ 6 , 
are consistent with Standard Model predictions, and can be translated into rough 
(conservative) limits on B + — > K + a by requiring that this exotic decay not exceed 
the total measured rate (we focus on the decays to K rather than K* because the 
observed rates are lower). The BaBar measurement includes lower-mass electron pairs, 
down to 30 MeV (compared to 140 MeV at Belle). The most recent Belle analysis 
[54] bins events by invariant mass, so that we can obtain a tighter limit (« 10~ 7 ) on 
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B + — > K + a in the region of interest, 140-1440 MeV. In all cases we take the limit to 
be the central value plus 2a. The BaBar measurement required that the £ + £~ pair 
originate from the same vertex as the K + . To set a conservative limit we require 
that the a decay within 100/mi, the scale of BaBar's vertex resolution. For the Belle 
analysis we require only that the PNGB decay within 0.5 cm, which would pass the 
requirement of [59J. 

• Similar but slightly weaker limits are obtained from measurements of K — > ir£ + £~, 
e.g. [60TI62] . We refer the reader to the original results and to [43] for details. 

It is worth emphasizing that the crude limits we have obtained are far weaker than the 
tightest limits that could be obtained by a directed analysis of BaBar or Belle data. Firstly, 
much tighter limits could be obtained by binning the 20-100 events in each sample more 
finely, and accounting for the detection efficiency as a function of mass. Further improvement 
could be obtained by including more displaced a decays, at the edge of the inner tracker (3 
cm) or even beyond. This direction is particularly worthy of exploration for the B — > Kfi + fi~ 
mode, for which there are no complementary searches near the high-F boundary of the Belle- 
excluded region, and for which the muon system gives an additional handle for studying highly 
displaced decays. 

B. Limits from Supernova SN 1987a 

For completeness we also include the constraints on PNGBs from SN 1987a [63]. We adapt 
the analysis of |64j to our setup, and obtain limits based on the assumption that PNGBs 
must not be the dominant mechanism of energy loss from the supernova. The temperature 
of the supernova core is conservatively estimated at T ~ 30 MeV [M], so PNGBs with 
mass significantly greater than these energies cannot be produced by the supernova, and are 
therefore unconstrained. The flux of PNGBs from the core of the supernova is approximately 



However, if these PNGBs decay or are re-absorbed, then they will not escape from the 
supernova and so they will not carry away any energy. Scattering and re-absorption dominate 



[64j 
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FIG. 2: Left: Constraints on pseudo-Nambu-Goldstone bosons as a function of their decay 
constants F and their mass m a from various meson decays: K + — >■ anything + e + e~ (green), 
-»• 7T++ invisible (blue), 5+ -»• K + £ + £~ (yellow) (£ = e,/x), and B + -> K + + invisible (red). 
Constraints from T(1S) or T(3S) — > — > 7+ invisible and K + — > ir + £ + £~ decays are weaker 
than those from B + — > K + + invisible and B + —> K + £ + £~, respectively, and thus not shown. 
Details are in §11 A Right: Gray shaded background region is the combined exclusion region from 
meson decays in the left figure. In the green exclusion region, the proton beam dump experiment 
CHARM at CERN would have seen at least five events (this exclusion region agrees roughly with 
that in [2]) - see [TV Here the PNGB is produced directly in the proton dump by a small mixing 
with the pion. For m a < 2m^, the PNGB decays to an electron pair, while in the "bubble" for 
m a > 2m.n the PNGB decays predominantly to a muon pair. The blue region is the limit from the 



supernova SN 1987a (see ^IIB). The light red region is the constraint from the muon anomalous 



magnetic moment and fills the gap for low m a and F left by the meson constraints (see ^IIC). The 
region excluded by the Fermilab E137 dump lies mostly within the CHARM excluded region and 
is not shown (it is instead shown in Fig. ^J. 



over PNGB decay for the relevant region of parameter space, giving a mean free path of 

A *~ 10m (io4v) 2 (8) 

for PNGBs with F < 10 8 GeV. We see that for F significantly smaller than 10 6 GeV, 
the PNGB mean free path is much less than the estimated core size of 10 km, so for these 
smaller values of F, SN 1987a does not constrain the PNGB. The exclusion contours we have 
plotted in Fig. [2] (right) correspond to requiring that the PNGBs carry away less energy than 
neutrinos, meaning that the total integrated PNGB emission must be less than about 10 53 
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PNGBs, each with energy of order T. 



C. Limits from the anomalous muon magnetic moment 

PNGB's contribute to the anomalous magnetic moment of the muon, a M , at the loop 
level. For the mass range of interest in this paper (< 1 GeV), only the one-loop contribution 
is important, and it is given by (see e.g. 



i K f 1 , x 3 



< = n — r / dx (9) 

8tT 2 F 2 J X 2 + (1 _ X )4 1 > 

\ / raj, 

M 

While the experimental measurement of is rather precise [66], the Standard Model pre- 
diction involves a hadronic contribution that must be estimated from experiments, which 
do not all agree. Using data from e + e~ annihilation to hadrons, the theoretical value of 
a M is smaller than the measured value by (316 ± 79)xl0~ n [67J, which is a 4cr discrepancy. 
However, estimates from r's give a smaller disagreement, with [68J finding a difference of 
(157 ± 82)xl0~ n , which is a 1.9cr discrepancy. 

Since the contribution from PNGB's is negative, a very conservative limit is obtained by 
using the 5a lower bound in [68], i.e. 

< > (157- 5 x 82) x 10" 11 , i.e. 

a£ > -253 x 10" 11 . (10) 
This constraint is included in Fig. [2] (right) and [5] 



III. CONSTRAINTS FROM LSND ON PNGBS AND DARK GAUGE BOSONS 

The Liquid Scintillator Neutrino Detector (LSND) experiment ran at the Los Alamos 
Neutron Scattering Center (LANSCE) in the 1990's [69J, and dumped C(10 23 ) 800 MeV 
protons on a predominantly water-copper target. This produced ~ 10 22 pions, a very large 
number that allows LSND to be sensitive in principle to very weakly coupled PNGBs or 
gauge bosons. A detector of length 8.3 m and a diameter of 5.7 m was located 29.7 m away 
from the target, 12° off-axis [84] , and filled with dilute liquid scintillator (there was no open 
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-Ebcam (GeV) 




X t (m) X d (m) 




E a (GeV) 


CHARM [2] 


400 


2.4 x 10 18 


480 515 


0.12 


25 


LSND (69H7T] 


0.8 


~ 10 23 


29.7 38 


see text 


0.3 


MINOS / MINERvA [TJ E 


g 120 


3.8 x 10 20 


1050 1087 


0.0006 


20 


MiniBooNE M 


8.9 


10 21 


541 553 


0.002 


2.7 



TABLE III: Shown are the total number of incident protons iV„, the distance from the target to 
the open decay region in front of the detector (i.e. the thickness of the shield) Xt, the distance from 
the target to the end of the detector Xd, the geometric acceptance times the number of pions per 
incident proton n^oegeo, and the median BNGB energy E a . These numbers were used to calculate 
the sensitivity of CHARM, MINOS /MINERvA (we always use the larger MINOS detector for 
estimates), and MiniBooNE to PNGBs produced directly in the proton dump. For LSND, we used 
a more involved procedure described in the text. 



decay region in front of the detector). Although the LSND collaboration did not search 
for signals that could originate from decays of long-lived exotics, approximate limits can be 
extracted from two published LSND analyses. 

We begin by reviewing the production of PNGB's in proton-nucleus collisions (which 



will re-appear in our discussion of more recent neutrino experiments in [TV), then discuss 
the implications of two specific LSND analyses [TOJ EI] f° r PNGBs that decay to e + e~ 
pairs or [i + pairs, respectively. The calculation involves considerable uncertainties and 
assumptions, and should be taken only as a sensitivity estimate — a dedicated analysis 
by the LSND collaboration is required to obtain a reliable limit. However, this analysis 
demonstrates that a dedicated analysis by the LSND collaboration would set some of the 
tightest constraints on PNGBs that couple to both quarks and leptons. 

We next discuss the production of dark gauge bosons in pion decay, and the implications 
of the same two analyses for these models. Our results are consistent with [7J, but we have 
considerably elaborated the discussion of the experimental sensitivity. The LSND sensitivity 
for this model overlaps closely with that of the SLAC electron beam-dump experiment E137 
0], as computed in [6j. 



A. Production of PNGBs at LSND and other Proton Beams 

In this section, we will consider the production and experimental sensitivity to pseu- 
doscalars from proton beam dumps. We will focus on the LSND experiment because it can 
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set the most stringent limits, although the MINOS and MiniBooNE experiments can also 



set interesting limits. We will discuss them more extensively in [TV 

Proton beam dumps can produce pseudoscalars directly through their mixing with pions. 
If the pseudoscalars couple to quarks they will interact via the operator 

— ±aqq =^ c— ^avr u (11) 

where c is an 0(1) parameter that depends on the up and down quark masses and any 
coefficients in the pseudoscalar coupling to the quarks. Since the pseudoscalar mixes with 
the pion, for every pion that is produced through a QCD process there is a probability of 
approximately (for c = 1) 

(t)' < 12 > 

of instead producing a PNGB. 

We can estimate the production rate within the detector acceptance and the PNGB 
momentum distribution by using the measured rates for it : 

F n 2 



N a = yyrj n n oN p e geo . (13) 

Here, n n o is the number of pions produced per incident proton, N p is the total number of 
protons, and e geo is the geometric acceptance (the solid angle subtended by the detector at 
the target divided by the solid angle of the beam). In practice, only the product n n o N p e geo 



is relevant but we have tabulated estimates for n w o e geo and N p in Table III for various 



experiments, for the sake of comparison. To determine the number of observable e + e~ pairs, 
we must also account for the probability of decaying in the detectable region (either inside 
the detector or in an open region upstream of the detector): 

N e = N a I — e~i^ ) ~ iV a e - ^, (14) 

where X t and Xd are the minimum and maximum decay lengths (roughly, X t is the thickness 
of the shield and Xd is the distance from the target to the end of the detector), and in the 
second expression we have assumed that jct ^> Xd — X t . 
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10" 2 0.1 1 

m a (GeV) 

FIG. 3: Sensitivity of various neutrino experiments to pseudo-Nambu-Goldstone bosons as a func- 
tion of their decay constants F and their mass m a . The thick (thin) black solid line corresponds 
to 10 (1000) events in LSND, the thick (thin) dashed blue line corresponds to 3 (1000) events 
in MiniBooNE, and the thick (thin) dot-dashed red line corresponds to 3 (1000) events in MI- 
NOS/MINERvA (in each case the inner regions correspond to more events than indicated by the 
line). Here the PNGB is produced directly in the proton dump by a small mixing with the pion. 
For m a < 2m the PNGB decays to an electron pair, while in the "bubbles" for m a > 2m „ 
the PNGB decays predominantly to a muon pair. The gray shaded regions are the combined ex- 
isting constraints from other beam dump experiments, meson decays, anomalous muon magnetic 
moment, and SN 1987a shown in the right plot of Fig. [2] 

B. LSND Analyses Sensitive to PNGBs 

We now focus on two LSND results and their implications for PNGB's. 

PNGBs decaying to e + e~ : The analysis in [70j used ~ 0.92 x 10 23 protons on target and 
looked for z/ M — > v e oscillations using from 7r + decay in flight [85]. The v e are detected 
via the inclusive charged-current reaction v e + 12 C — > e~ + X. This analysis focused on 
identifying electrons in the energy range 60 MeV to 200 MeV. Various cuts were used in the 
analysis with an energy-dependent efficiency that is always near 10%. Clearly this analysis 
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should be sensitive to PNGBs decaying to electrons inside the LSND detector, although it is 
impossible to accurately estimate the efficiency without a dedicated analysis. One difference 
between our signal and the study in [70] is that in our case both an electron and a positron 
are produced in the detector, as opposed to just a single electron. However, it is impossible 
to distinguish e + e~ events from a single electron (or single photon) event, and it has been 
suggested [75] that we should assume that the total e + e~ pair energy (i.e. the PNGB energy) 
would have been measured as the energy of a "single-electron" . For simplicity, we will assume 
that the detection efficiency for an e + e~ pair is roughly the same as for the single electron 
analysis, i.e. e^i ~ 0.1. 

We have to estimate the energy distribution and number of PNGBs incident on the LSND 
detector. This should be roughly equivalent to the number and energy distribution of pions, 
which has been simulated with Monte Carlo by the LSND collaboration [75]. Specifically, 
we will model the PNGB kinetic energy distribution using the predicted pion kinetic energy 
distribution from [7H] , and rescaling the total rate by F 2 / F 2 from the total number of pions 
incident on the detector, (8.6 ±2.1) x 10 14 cm -2 . This should be an excellent approximation 
for m a « m n ; it may be subject to additional 0(1) uncertainties for m a <C and especially 
for m a ^> m n , which we neglect. In addition to the reconstruction efficiency, we must also 
account for the fraction of PNGBs with kinetic energy between 60 — m a and 200 — m a , which 
for small m a and the assumed distribution is approximately 25%. We note that the mean 
of the pion kinetic energy distribution is at ~ 275 MeV with a root-mean- square spread 
of ~ 130 MeV, so an analysis including higher-energy electrons would be significantly more 
sensitive. 

We show the number of e + e~ events obtained from PNGB decays in the F versus m a 
parameter space in Fig. [3] where the solid black thick and thin lines for m a < 2m M show 
10 and 1000 e + e~ events, respectively. We have assumed that the number of PNGB decays 



inside the detector is given by Eq. ( 14 ) integrated over an PNGB energy from 60 MeV — m a 
to 200 MeV — m a , and multiplied by the efficiency e e s t i- The analysis of [70] (see their 
Fig. 29) indicates that 10 PNGB events in a 20 MeV energy bin below 200 MeV would have 
been easily noticed. We note that the sensitivity could have been increased by increasing 
the energy threshold for electrons beyond 200 MeV. 

PNGBs decaying to /i + yT : We now turn our attention to the analysis in [7T], which would 
have been sensitive to events. The analysis in [7T] considered ~ 1.8 x 10 23 protons on 
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target and searched for 7r° — > u^u^. The neutrinos interact with the nuclei in the detector 
to produce muons through the reactions + 12 C — > yT +p + X and + 12 C — > fi + + n + X. 
This analysis focused on identifying high-energy, muon-like beam excess events in the energy 
range 160 MeV to 600 MeV electron equivalent. The muon was required to decay inside 
the detector. Various selection cuts were used in the analysis [71] with an overall efficiency 
of ~ 0.15 for identifying neutrino interactions in the detector. This analysis should clearly 
have been sensitive to PNGBs decaying to yU + /i~ pairs inside the LSND detector, although 
it is impossible to accurately calculate the efficiency without a dedicated analysis. We 
simply estimate the efficiency for n + [i~ pairs to be reconstructed as single-muon events to 
be e c fi t 2 ~ 0.1, a little bit less than the efficiency for the original analysis. We estimate the 
PNGB kinetic energy distribution as described above. To obtain the total number of pions 
incident on the detector, we rescale the value we used above for the a — > e + e~ analysis to 
1.7 x 10 15 cm -2 , since this analysis uses a factor of ~ 2 more protons. Integrating over an 
PNGB energy from 160 MeV — m a + 2m^ to 600 MeV — m a + 2m M , we show the number 
of events obtained from PNGB decays in the F versus m a parameter space in Fig. |3j 

where the solid black thick and thin lines for the region m a > 2m M show 10 and 1000 fi + fi~ 
events, respectively. The analysis in [71] indicates that 10 PNGB events in a given energy 
bin would have been noticed, and so we believe that this is a reasonable number of events 
with which to set a tentative limit. 

C. LSND Limits on Light Vector Bosons 

In [7J it was shown that LSND may be sensitive to a new light vector boson A' that 
mixes with U(l) hypercharge; however, the precise sensitivity was unclear due to uncertain- 
ties in the LSND experiment. We believe that we have clarified these uncertainties using 
information from [701 [7TJ [75] . 

With approximately 1.8 x 10 23 protons on target at 800 MeV, the LSND experiment 
produced about 1.4 x 10 22 7r°'s. The vast majority of these 7r°'s decay to 77, so with a 
kinetic mixing parameter e we expect a branching fraction of 2e 2 for 7r° — > ^A'. We will only 
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consider the simplest scenario, where the A' decays to e + e with a lifetime 

that is also set by the kinetic mixing. Due to the huge number of pions produced at LSND, 
the A' may be visible even for very small e. 

To set a reliable limit, we use the results of the n° — > vv search discussed previously [7Tj . 
They performed a Monte Carlo simulation to determine the flux of neutrinos in the LSND 
detector from this decay and plotted it in their Fig. 1 [71]. In our case we want to consider 
7T° — > ^A', but in the limit that <C the neutrino results should give a very good 
approximation to the distribution of A' particles; beyond this limit we will continue to use 
the neutrino results with the caveat that we expect an 0(1) uncertainty in the rates. Given 
this distribution of A' particles, we computed the number of A ; s that would have decayed 
inside the LSND detector using the lifetime formula above, and the 2e 2 branching ratio of 

The only reliable data that can exclude an A' decaying to e + e~ inside the LSND detector 
is the analysis of |70| . which we discussed above in order to determine the LSND sensitivity 
to PNGBs. From the study of tt° backgrounds in [70] and from discussions with LSND 
collaborators [75], we expect that an A' decay to e + e~ would have been interpreted as a 
single-electron event in the LSND detector. The data on these types of events has only been 
published for energy depositions below 200 MeV, and the analysis in question had an overall 
reconstruction efficiency of about 10%. Thus we use the Monte Carlo data [TQl [75] below 
this energy to set our tentative limits, which are displayed in Fig. |4} Our limits are based 
on the conservative assumption that a signal of 50 or more events after the 10% efficiency 
cut should have been visible. 



IV. SENSITIVITY OF MODERN NEUTRINO EXPERIMENTS 



A. Production from Proton-Nucleus Collisions 



We have discussed this production mode already in the context of LSND, in §111 A A 



similar analysis was also carried out by the CHARM experiment at CERN [2], a proton 
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FIG. 4: The sensitivity of various beam dump, collider, and astrophysical probes of light vector 
bosons that kinetically mix with the standard model hypercharge, as a function of the kinetic 
mixing parameter e and the vector boson mass rriA' (from [6]). The thick (thin) solid black line 
corresponds to 10 (1000) signal events in LSND. The region enclosed by the thick line can be viewed 
as a very rough exclusion limit from the LSND experiment. A re-analysis of the LSND data by the 
LSND collaboration could further extend the sensitivity of that experiment [7J. Further details are 
described in the text. 

beam dump with a 400 GeV beam, a 35 m decay region, and a 3 x 3 meter detector. Similar 
searches could be done using the MINOS or MINERvA detectors on the NuMI beamline, 
and with the MiniBooNE detector. We summarize the parameters for these experiments in 
Table EH 

In all cases we estimate the production of PGNBs by rescaling the pion production rate 
(determined either experimentally or from Monte Carlo) by a factor of (F^/F) 2 . However, 
the production rate of very weakly coupled PGNBs could actually be significantly greater 
than this scaled rate because PGNBs produced within the target will virtually always escape 
from the target, whereas many of the pions produced by proton interactions might not make 
it out of the beam dump. This issue may merit further exploration. 

For the CHARM experiment, we take the experimental parameters from [2]. Our calcu- 
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lation of the PNGB exclusion region is shown in Fig. [2] (right) and roughly agrees with the 
results in [2J. 

The MiniBooNE estimate is shown as a thick (thin) blue dashed line in Fig. [3j corre- 
sponding to 3 (1000) events. It is obtained as follows. We take the 7r° production cross 
section to be the average of the ir + and 7r~ production cross sections found in [74J. These 
cross sections are given as a function of angle and momentum, and one can show that the 
fraction of 7r°'s produced in the target that point towards the 12 m detector at a distance 
of 541 m away is roughly 0.002. In addition, [73] gives the average number of tt + and 7r~ 
produced per "particle-producing reaction", which we take to be per incident proton on 
target - we average these to find about 0.89 7r 0, s per incident proton. 

The MINOS/MINERvA estimate is shown as a thick (thin) red dot-dashed line in Fig. |3j 
again corresponding to 3 (1000) events. This estimate for MINOS was obtained by using data 
on their 7r + momentum distributions from Fig. 9.3 of [72]. Based on the assumption that the 
momentum distribution for a general PGNB would be similar, we estimate from that figure 
that a fraction 0.0035 of PGNBs produced by the NuMI proton dump would point towards 
the MINOS detector. Data from [72] shows that there were roughly 0.18 pions produced per 
proton on target, which allows us to estimate the total number of PGNBs produced. 

Searches in either MiniBooNE and MINOS/MINERvA would have significant overlap 
with the existing CHARM result, perhaps slightly extending the region probed. However, 
the proposed "Project X" upgrades to NuMI, which would increase the beam intensity by 
a factor of ~ 5-10 and possibly include a larger near detector could extend sensitivity to 
larger F. 

B. Production through Muon a-sstrahlung 

A more unique search opportunity at MINOS/MINERvA and MiniBooNE arises from 
their magnetic focusing of charged pions, which in turn focuses the neutrinos and muons 
from their decay towards the detectors. As a result of this focusing, approximately one in a 
thousand protons on target produces a muon that points toward the neutrino detector. As 
these muons stop in the rock upstream of the detector, they can radiate very forward PNGB's 
in a process that is analogous to ordinary bremsstrahlung. Unlike photon bremsstrahlung, 
however, the typical PNGB produced by this process carries a large fraction of the incident 
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£beam (GeV) Nj E e (GeV) X t (m) X d (m) W (m) E thresh ( GeV) 

E137 jl] 20.0 N e = 1.2 x 10 20 20 180 380 3 2.0 

MINOS/MINERvA [76] 120.0 = 2.7 x 10 17 20 240 270 3 1.0 

MiniBooNE [77] 8.9 = 2.0 x 10 17 1.3 450 462 12 0.1 

TABLE IV: Shown are the proton beam energy -Ebeamj the total number of incident electrons (for 
E137) or muons (for MINOS /MINERvA and MiniBooNE) N e and their average energy E e , the 
distance from the start of the electron or muon dump to the open decay region (if any) in front of 
the detector or to the detector itself (i.e. the thickness of the shield) Xt, the distance from the start 
of the muon dump to the end of the detector X^, the diameter of the detector W, and the threshold 
energy -Ethresh to detect an electron or muon that originates from an PNGB decay. These numbers 
were used to calculate the sensitivity of E137 |4|, MINOS/MINERvA {76] (we always use the larger 
MINOS detector for estimates), and MiniBooNE [77] to PNGBs produced by bremsstrahlung off 
an electron beam (in the case of E137) or off a muon beam that is produced in a proton dump. 



muon's energy. The production rate and kinematics can be reliably calculated using the 
Weizsacker- Williams approximation [75], where the nuclei in the fixed target provide an 
effective photon beam. We have relegated the details to Appendix [A] so in this section we 
will simply give the results along with their physical motivation. 

We will be interested in a wide range of values for m a ; with this in mind we note that 
production is dominated by emission angles 

a < max ( !£, ^i) (16) 



for E a ~ Eq. This is very useful for estimating the angular acceptance, which we compute 
by taking the ratio of the solid angle subtended by the detector to the solid angle within 
6 a of the beam direction. We also considered acceptance issues associated with the PNGB 
decay, but because we require the PNGB to decay within the MiniBooNE detector or in the 
relatively small open decay region in front of the MINOS/MINERvA detector, this is not 
important, and we ignore it. 

The cross section is always peaked near E a pa Eq, even when m a «C m^. This is in stark 
contrast to the analogous formula for photon bremsstrahlung, where the rate is proportional 
to the inverse of the photon energy. The difference is due to the contrasting soft emission 
behavior of gauge bosons and goldstone bosons - the former have soft singularities, while 
emission of the latter vanishes in the soft limit due to their derivative couplings. The 
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FIG. 5: Sensitivity and constraints of various experiments to leptophilic pseudo-Nambu-Goldstone 
bosons as a function of their decay constants F and their mass m a . Here the PNGB is produced by 
bremsstrahlung off an incident muon or electron beam. Thick (thin) lines show rough sensitivity 
regions and correspond to 3 (1000) displaced e + e~ pairs in MINOS/MINERvA (red dot-dashed 
lines), MiniBooNE (blue dashed lines), and in a thin target experiment using the COMPASS muon 
beam (green dotted line, see j]v|). The thick (thin) dotted black lines correspond to S/\ /r B = 3 (10), 
where S (B) are the number of prompt (jP~fJ,~ signal (radiative background) events in COMPASS 
(see § |v| ) (we have ignored the Bethe-Heitler background and the finite acceptance, so these lines 
should not be viewed as real significance lines but only as very rough estimates of what could be 
probed). Inside the gray shaded region, E137 would have seen at least one event - since they saw 
none, this region gives their approximate constraint. Details are described in the text. The light 



red region is the constraint from the muon anomalous magnetic moment (see { II C ) . 



peaking toward high energy fraction E a /Eo is further enhanced by a larger phase space 
when m a ^> m^. The total cross section for pseudoscalar production from a muon beam has 
the parametric form 

t A max(m^, m A a ) 

Note that the formula in the case of an electron beam dump would be identical except with 
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m M — > m e , so we see that muons are a much more efficient source of PNGBs as long as 
m a > m e . 

To get a rough idea of our experimental reach, it is useful to have approximate formulas 
for the total production rates. For a thin target, the yield is 

N a ~N< T f$ (18) 
at 1 max(m~ , m^j 

where T e is the number of electron radiation lengths of material. In the case of a thick 
target, where the beam of muons is completely dumped, it is more difficult to give a simple 
parametric formula for the pseudoscalar yield because muons stop due to minimum ionization 
interactions, as opposed to bremsstrahlung. However the equation above with T e ~ 100 
provides a rough estimate of the pseudoscalar yield for MINOS and MiniBooNE. 

Fig. [5] shows the number of e + e~ pairs in MINOS/MINERvA (red dot-dashed lines) and 
MiniBooNE (blue dashed lines) as a function of the PNGB decay constant F and the PNGB 
mass m a . The thick (thin) lines correspond to 3 events (1000 events). The gray shaded region 
corresponds to the approximate constraint from E137: inside the region, more than 1 e + e~ 
event would have been seen (note that we calculated this region using the procedure in [6], 
but changing the couplings from an A' to those relevant for PNGBs). Our estimate for 
the E137 constrained region agrees with that in jl]. We see that MINOS/MINERvA can 
extend the E137 region, although the MiniBooNE region is contained within the E137 region 
since many PNGBs are produced in the dump with a large enough angle causing them to 
miss the detector (this is because the MiniBooNE muon beam has a lower energy and their 
detector is further away compared to MINOS/MINERvA). Note that these lines have been 
calculated with the full formula as detailed in the appendix. The NOVA experiment and 
"Project X" upgrades to the NuMI beamline will have a factor of 5-10 more protons and 
so will significantly extend the reach. The light red region is the constraint from the muon 
anomalous magnetic moment, while the other lines will be discussed in §|V) 
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V. MUON FIXED-TARGET EXPERIMENTS WITH THIN TARGETS COM- 
PASS 



In this note, we have argued that new light states may be produced by bremsstrahlung 
off a dumped muon beam. Here we will briefly comment on the potential of using a fixed 
target setup to constrain weakly coupled light states. Compared to a dump experiment, 
the target will be much thinner, and we can thus search for particles with a much shorter 
lifetime. 

Electron fixed target experiments have been used to probe new light states |JJ S3 EH] , 
but muon fixed-target setups are different in several respects. The intensity of muon beams 
is obviously lower, but this may be partly compensated for by a much thicker target - 
a muon beam can easily traverse a meter of material leaving a relatively quiet off-beam 
environment. Furthermore, muon beams will have an obvious advantage over electrons in 
PNGB searches, as their couplings are proportional to particle mass. For example, it is 
instructive to consider the COMPASS experiment at CERN [80], where a 160 GeV muon 
beam strikes a low-Z polarized target (they have also used higher-Z, but much thinner, 
targets). A detailed description of their polarized target is given in [80], and we approximate 
it as a 130 cm long target consisting of Lithium, with a packing factor of 0.5 (i.e. a column 
density of 34.7 g/cm 2 , or 0.42 radiation lengths). In total, we estimate that this experiment 
has collided about ~ 10 15 muons. Here we will illustrate the rough reach of this COMPASS- 
like setup assuming the Lithium target. We have checked that a higher Z target, such as 
Tungsten, of a similar thickness would probe a somewhat larger parameter space, including 
high-F regions above the muon threshold. We will not show this in any figures, focusing 
instead on the data set already collected by COMPASS. 

There are two regions of parameter space to consider, one in which F is low and the PNGB 
decays promptly, and one in which F is high and the PNGB decays with a displaced vertex. 
In the case of a prompt decay, one must search for a peak in the di-muon invariant mass on 
top of a sizable standard model background (decays to electrons may not be searched for 
this way because their interactions in the target considerably degrade the mass resolution). 
The signal to background ratio (where the background comes from di-lepton production via 
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an off-shell photon) for such a search is roughly [6] 

S ml 3 m a 

— ~ — (19) 

B F 2 Aa 2 5m' 1 ; 

where 5m is the resolution-limited mass window used in the search. Based on the angular 
and momentum resolutions for COMPASS [81], we infer a mass resolution a m = 11 MeV, 
making a window width 5m ~ 2.5<r m = 27 MeV appropriate for estimating sensitivity. 
Additional two-photon diagrams (which we refer to as the "Bethe-Heitler" background) and 
their interference have been neglected; these actually dominate over the "radiative" dilepton 



production assumed in [T9j but can be reduced to an order-1 fraction of the background by 
kinematic selection. 



We can combine Eq. (19) with the expected number of signal events using the formulas 
in Appendix [A] to get a very rough estimate for the sensitivity of such a search. These 
estimates omit several factors, all of which degrade sensitivity: the finite detector acceptance, 
the additional Bethe-Heitler background, and the finite acceptafinite of kinematic selection 
needed to veto the dominant Bethe-Heitler background. Assuming 10 15 muons on target 
and a mass window 5m = 27 MeV [81], we show the "3<t" and "10<r" estimated reach in 
Fig. [5] (the neglected factors are expected to reduce the sensitivity each point in parameter 
space by a factor of 2-4). Such a search in existing COMPASS data (or that of a similar 
experiment) would be sensitive to leptophilic PNGBs with F near the weak scale and masses 
between the dimuon threshold and a few GeV. Higher F's can be probed by increasing the 
target thickness in radiation lengths, for example by using a high-Z target of comparable 
length. 

Another parameter range that can be studied by muon fixed-target experiments is that 
where the PNGB decays are significantly displaced. If the vertex can be reconstructed to 
be downstream of the target region, beyond the tails of Standard Model backgrounds, a few 
events could be enough to claim a discovery. The transverse vertex resolution of COMPASS 
is about 0.1 mm [80], so that the vertex resolution in the direction of the beam is about 
o z = (0.1 mm) x E heam /m a . In Fig. [5j we show the area of parameter space that yields 3 
or 1000 events in the region between 5a z and 10 meters behind the 130 cm Lithium target 
(for this estimate we assume that the PNGB is produced in the middle of the target). Note 
that the vertex resolution is worse for lower-mass PNGBs, so that for m a < lOMeV, we 
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have 5er z > 10 m, and there is no good search region (which is why the plot cuts off for low 
m a ). An analysis of the data in COMPASS, or a similar experiment of this type, can cover 
new ground with respect to E137 or the neutrino experiments at lower F, and can close the 
window between the latter searches and the muon g — 2 limit. With a long higher- Z target 
such as Tungsten, more PNGBs would be produced so that some region of parameter space 
at F ~ 10 4 GeV above the muon threshold may also give rise to observably large rates of 
displaced decays. 

In summary, the very rough estimates above suggest that a fixed target experiment with a 
focused muon beam may be able to probe unconstrained regions in PNGB parameter space. 
In fact, the existing data set of the COMPASS experiment may already be able to set some 
interesting new limits for leptophilic PNGBs. It may also be worthwhile to consider using 
more diffuse muon "beams" such as those in neutrino factories in a thin target setup due to 
their higher intensity. We leave this for future thought. 

VI. CONCLUSIONS 

We have explored the sensitivity of neutrino oscillation experiments to three types of new 
light states - vector bosons that kinetically mix with the photon, pseudoscalars that couple 
to quarks and leptons, and pseudoscalars that couple preferentially to leptons. The first two 
are strongly constrained by rare decays, fixed target experiments, and supernova 1987a (all 
of which we have reviewed), whereas there are fewer tests of the third class. 

The sensitivity of the LSND experiment to vector bosons was discussed in [7], but the 
details of the LSND detector and analyses were not considered. We have shown that the 
analyses of [701 EI] would have been sensitive to vector bosons with mass below 2m M and 
a large range of coupling strengths. This sensitivity has significant overlap with the E137 
experiment [I] , but LSND does probe a new region at very weak coupling, and in any case the 
LSND results serve as an important cross-check. Because LSND dumped a larger number of 
protons, other neutrino experiments are not as sensitive to the light vector boson scenario. 

We also considered the sensitivity of various neutrino experiments to pseudo-scalars that 
couple to quarks, so that they can be produced in proton beam dumps. LSND remains the 
most sensitive experiment for pseudoscalars with mass below 2m M , nearly closing the gap 
between fixed target experiments and supernova constraints. However, other experiments 
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are more sensitive to heavier pseudoscalars. The MiniBooNE and MINOS experiments are 
currently competitive with the best limits on these particles. However, the estimates that 
we derive for these experiments may be too conservative because we estimate the rate of 
PGNB production by scaling the pion production rate, which may be an underestimate if 
many pions get stuck in the beam dump. This issue merits further study. 

The MINOS and MiniBooNE experiments produce neutrinos from focused muon beams; 
the requisite muon beam dumps provide a unique opportunity to search for pseudoscalars 
that couple preferentially to leptons, since we expect these particles to couple far more to 
muons than to electrons. The MINOS experiment is sensitive to leptophilic pseudoscalars 
with decay constants almost an order of magnitude greater than any previous experiment 
(in particular E137), while MiniBooNE can probe a region that is contained within that of 
E137. 

A thin target experiment with a muon beam, such as that available in COMPASS, offers 
a unique probe for leptophilic PNGBs. An analysis using the existing COMPASS data set 
and looking for either e + e~ pairs originating from displaced vertices behind their (Lithium) 
target or for a spike in the fi + fi~ invariant mass spectrum of muons coming from the target 
should be sensitive to new regions of parameter space. A similar experiment using a higher- Z 
target would have even more sensitivity. 

Upgrades to the NuMi beamline followed by the proposed "Project X" experiment will 
explore new parameter space for both standard and leptophilic pseudoscalars. It is our hope 
that in the future these experiments will perform dedicated analyses to explore and constrain 
new weakly coupled low-mass particles. 
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Appendix A: Pseudoscalar Production 



We are interested in the pseudoscalar production rate from muons braking in a fixed 
target. This process is analogous to ordinary bremsstrahlung, and it can be reliably calcu- 
lated using the Weizsacker- Williams approximation [78], where the nuclei in the fixed target 
provide an effective photon beam. When the incoming muon has energy Eq, the differential 
cross section to produce a pseudoscalar of mass m a with energy E a = xE is 



da 



dx dcos9 n 
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(A2) 



is the virtuality of the intermediate muon in initial state bremsstrahlung, and x is a form 
factor that can be found in [78J. We will be interested in a wide range of values for m a ; with 
this in mind we see that production is dominated by 



. m a m M 
a < max ( — , — 
tin £jq 



(A3) 



for x ~ 1. For angles larger than this, the differential cross section falls off rapidly, as 1/0*, so 
this angular scale sets the width of the PNGB beam for the purposes of angular acceptance. 
We compute the angular acceptance by taking the ratio of the solid angle subtended by the 
detector to the solid angle within a of the beam direction. 



Integrating Eq. (Al) over the angle a , we find 
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and we have now included the form factors 



Xi(Z) = Z 2 ln(184Z" 1/3 ) + Zln(1194Z~ 2/3 ) (A6) 
X 2(Z) = Z 2 + Z. (A7) 

It is important to note that the cross section is always peaked near x ~ 1, even when 
m a <C m M . This is in stark contrast to the analogous formula for bremsstrahlung, which is 
proportional to 1/x. The difference is due to the different soft emission behavior of gauge 
bosons and goldstone bosons - the former have soft singularities, while emission of the latter 
vanishes in the soft limit. The cross section is dominated for x ~ 1 with the parametric 
form 

m l 2a 2 

a -~d — r~2 — A8 

t A max(m^, m A a ) 

The max(m 2 ,m 2 ) factor comes from the presence of the function / ~ m 2 /m 2 in the de- 
nominator for larger m a /m^. The formula in the case of electrons would be identical with 
m M — > m e , so we see that as claimed, muons are a much more efficient source of PNGBs as 
long as m a ^> m e . 

To use these formulae we must account for the way that the muon slows in a beam dump. 
The number of PNGBs produced per incident muon is 



^^f\ El f\ tM E aM % (A9) 



where E is the energy of the original incident muons, E a is the energy of the produced 
PNGBs, x' = E a /Ei, x = E a /E , Na is Avogadro's number, X e is the unit (electron) 
radiation length in g/cm 2 , A is the atomic number of the material, and T e is the total number 
of (electron) radiation lengths in the target or beam dump. The function I^Eq, E±, t e ) is 
the distribution of muon energies after the muons have traversed t e radiation lengths. 

It is essential to remember that muons stop primarily through minimum ionizing interac- 
tions, and not through radiation, so the number of radiation lengths is not directly related 
to the muon energy loss. We therefore use electron radiation lengths in our computations, 
so that t e and T e can be much greater than 1 without completely depleting the energy of 
the muon beam. We estimate the function 1^ using the relevant material properties (if the 
intervening material is earth, we use Silicon). 
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To compute the number of electron or muon pairs from PNGB decays in a detector, we 
simply integrate Eq. (A9) times the probability that an PNGB with this energy decays in 
(or in some cases in front of) the relevant detector (see Eq. (14)). 
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